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Abstract 

We investigate nonperturbative infrared effects for the 0{N) linear sigma model in de Sitter 

CO ' space using the two-particle irreducible effective action at the Hartree truncation level. This 

O . approximation resums the infinite series of so-called superdaisy diagrams. For the proper treatment 

Sh . of ultraviolet divergences, we first study the renormalization of this approximation on a general 

P-c 

■^ I curved background. Then, we calculate radiatively corrected masses and the effective potential. 

*vj ' As a result, spontaneous symmetry breaking is possible, on the other hand, the Goldstone modes 

acquire a positive definite mass term due to the screening effects of interaction. Possible infrared 

. . divergence is self-regulated by the mass term. Furthermore, there is a symmetry restoring phase 

^ I transition as a function of the Hubble parameter. In our approximation, the phase transition is of 

first order. 
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I. INTRODUCTION 

In de Sitter space, the propagator for a massless minimally coupled scalar field has an 
infrared divergence [l|. This causes a breakdown of perturbative expansion around massless 
fields. Moreover, in the theory with a sufficiently small mass, the propagator may have an 
indefinitely large term. In this case, perturbation theory can no longer be a valid approx- 
imation. Motivated by the inflationary cosmology, the interpretation and treatment of the 
infrared divergence has become of interest in recent years [jtJ]- Also, it suggests that the 
strong infrared effects may give novel effects on inflationary physics. 

One of the strategies for treating such light scalar fields is resummation of perturbative 
expansion. In this method, we include an infinite series of diagrams with certain topology. 
Then, it becomes possible to capture some nonperturbative quantum aspects, which cannot 
be obtained in ordinary perturbation theory. In this direction, the 0{N) linear sigma model 
is analyzed by using the large- iV expansion in the earlier study |5|. Furthermore, in our 
previous work, we study 0^ theory using the two-particle irreducible (2PI) effective action |6| . 
In both analyses, the fields acquire a positive definite mass term due to nonperturbative 
screening effects of interaction. The possible existence of an infrared divergence is self- 
regulated by the mass term. On the other hand, these analyses predict different behaviors 
for the effective potential. In the large-A^ analysis, it was shown that the strong infrared 
effects prevent the existence of spontaneous symmetry breaking. In contrast, for 0^ theory, 
the spontaneous symmetry breaking of Z2 symmetry is possible. Moreover, the broken 
Z2 symmetry is restored by the first order phase transition as a function of the Hubble 
parameter. 

The absence of symmetry broken states in de Sitter space is reminiscent of the result 
for scalar field theory in flat two dimensions. In fact, the propagator for a massless scalar 
field in flat two dimensions also has an infrared divergence. Then, the spontaneous breaking 
of continuous global symmetry is prohibited, avoiding the presence of massless Goldstone 
bosons |7|. This fact makes the absence of symmetry breaking in de Sitter space plausible. 
However, if it is true, there is a discontinuous transition between the broken phase in flat 
space and the symmetric phase in de Sitter space. Moreover, it remains open whether the 
result predicted by the large- A^ expansion is applicable to any number of A^. If a broken 
phase exists at some A^, the Goldstone modes appear to cause the infrared divergence. In 



such a circumstance, an investigation into whether the construction for well-defined theory 
requires some infrared cutoff or a natural regularization mechanism like the self-regulation 
exists is desirable. 

In this paper, we study the 0{N) linear sigma model by using the 2PI effective action 
at the Hartree truncation level. The 2PI effective action enables us to extract some non- 
perturbative quantum aspects by resumming perturbatively expanded vacuum diagrams |8| . 
However, it is a nontrivial task to renormalize ultraviolet divergences at a given truncation 
of the 2PI effective action. In fact, even at the lowest order truncation of the mean field 
approximation, renormalization is not trivial, where different renormalization prescriptions 
predict different results in the literature (9|, [lO| . However, recently much progress has been 



made regarding the renormalization of the 2PI effective action formalism. In Ref. ll|, a 
systematic treatment of ultraviolet divergences at a given truncation of the 2PI effective 
action for scalar field theory in fiat space is proposed. Furthermore, explicit counterterms 
are constructed for scalar field theory with more complicated symmetry at the 2PI Hartree 
approximation [12]. In our previous study, we derive the counterterms needed for 0^^ theory 



at the 2PI Hartree truncation on a general curved background [13|. In this paper, we extend 
this renormalization scheme to the 0{N) linear sigma model. Using this renormalization, 
we investigate the nonperturbative infrared effects on the mass and the effective potential. 

There is an earlier study with quite similar analysis, using the mean field approxima- 
tion 1^. However, we believe that the renormalization in this literature is not appropriate. 
This is because ultraviolet divergences are renormalized only by the mass counterterm, 
which depends on the vacuum expectation value and the Hubble parameter. Moreover, this 
inappropriate renormalization obstructs a further study in the framework of the Hartree 
approximation. A difference between our study and the preceding work is the renormaliza- 
tion prescription; our renormalization prescription enables us to study in detail at the 2PI 
Hartree approximation. 

This paper is organized as follows. In Sec. II, we provide detail as to how to renormalize 
the 0{N) model at the Hartree truncation level of the 2PI effective action. In Sec. HI, 
radiative corrections to the mass and the effective potential are calculated. Its similarity to 
finite temperature field theory is mentioned. Moreover, the consistency between our results 
and the large-A^ expansion is discussed. Sec. IV is devoted to a conclusion and discussion. 
In this paper, we use the unit system oi h = c= 1. 



II. RENORMALIZATION OF THE 2PI HARTREE APPROXIMATION FOR THE 
0{N) LINEAR SIGMA MODEL 

In this section, we show the renormahzation of the 2PI Hartree approximation for the 
0{N) model in general curved space. We first renormalize the theory in fiat space, then we 
present the way of extension to curved space. 

In flat space, the 2PI effective action for a multi-component scalar fleld, which is a 
functional of the vacuum expectation value of the quantum fleld v'' and the full propagator 



Gij, is given by 



a 



(1) 
where summation is understood in the repeated indices and 



r-^\v^](x X') - ^j31^^ 



iG'o-,M^1(x,x') = ---4-4-, (2) 



is an inverse propagator. 

r2[v^, Gij] is expressed by (— z) times all of the two-particle irreducible vacuum diagrams 
with the propagator given by G and vertices given by interaction terms S'int of the shifted 
action. (p{x) = 0(x) — v{x) is a shifted fleld. Here, a two-particle irreducible diagram 
is a diagram which cannot be cut in two by cutting only two internal lines, otherwise it is 
two-particle reducible. Various approximations can be made by truncating the diagrammatic 
expansion for T2[v, G]. The mean fleld and gap equations are given as a stationary condition 
for r[t;,G] with respect to v and G. From these equations, we can solve G as a function 
of f , G = G[v]. Then the standard IPI effective action is obtained by inserting G[v] into 
r[t;,G], giving ripiH = r[i;,GH]. 

For the 0{N) linear sigma model with the action 

5[0,] = -Jd'x[^<p\n + m2)0^ + ^(0>O'] , (3) 

where the index i runs from 1 to A^, the shifted action is given by the following 

S[ip^ + v'] = S[v'] + I d^xlg^'^d^^'d^v' - m^v'^' - -^4(t;V)(^;V^') - \\^ + ^^W 

J y 4iV 2 ^^^ 





FIG. 1. 2PI diagrams at the 2PI Hartree approximation. The wiggly hne represents the vacuum 
expectation value of the quantum field, f*. 

From the expression, we can read off the free inverse propagator and the interaction terms 
oi Ojjj^ as 



iGoi{x,x') 
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(5) 
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(/(/.*)^ + 4((/.V*)(^V') , (6) 



where we have defined GqI, to be independent of f*. Furthermore, the diagrams which are 
constructed from counterterm vertices with only one internal line are considered to be 2PI 
diagrams. 

We now approximate the theory by including only the tadpole and double bubble dia- 
grams as the 2PI diagrams as shown in Fig. [1] and the corresponding counterterm diagrams 
needed at this truncation. This truncation corresponds to the Hartree-Fock approximation. 

A. Renormalization 



In ordinary renormalized perturbation theory, counterterms are perturbatively expanded, 
and each counterterm vertices in the counterterm diagrams have the same expression at each 
order in an expansion parameter. However, in the resummed formalism, in general there are 
no expansion parameters. Then, each counterterm vertices in the counterterm diagrams do 
not necessarily have the same expression. That is, we have to introduce different counterterm 
vertices in each counterterm diagram. 

In our double-bubble approximation, the reason for introducing different counterterms 
is explained as the difference in the channels in the resummation process. At the two-loop 
order of the 2PI diagrams, there is the double-bubble diagram and the setting sun diagram 
as shown in Fig. 2. The double-bubble diagram has the contribution of the diagrams con- 



structed by the s-channel. On the other hand, the setting sun diagram has the contribution 
of the t- and w-channel diagrams. Inclusion of only the double bubble diagrams at the trun- 
cation means that we have discarded the contributions of the sub-diagrams constructed by 
the t- and u-channels. How the ultraviolet divergences of the discarded t- and u-channel 
diagrams have been engaged to counterterms is different for each counterterm vertex. Thus, 
we have to introduce different expressions to each counterterm vertex. 

In 0^ theory, we have to introduce different coupling counterterm vertices 6X for S[v] and 
T2[v, G]. In the 0{N) model, we have to further introduce independent counterterms to the 
counterterm diagrams of r2[w,G'] distinguished by the way of contraction. The reason is 
explained by the different role of each counterterm vertex diagrams in the renormalization 
process (see Appendix A). 

First, two tadpole diagrams as shown in Fig. [3] can be constructed from the two point 
interaction terms in Eq. ([6]). We assign the following independent counterterms for each 
contractions 



2iVL 



1 

2iV 



1 (7) 



where the superscripts, A and B, refer to the different counterterms as A"^ = A + ^A"^ and 
A^ = A + ^A'^, though the finite part is the same. 

From the four point interaction terms in Eq. (|6]), we can construct the two double bubble 
diagrams as shown in Fig. 3. Again, we assign the different counterterms as 



4JY [^0 2^ Gii 2_^ Gjj + 2A^ 2^ G 



(8) 



In the case of four point interaction, we assign the different counterterms by the different 
ways of contraction, though the interaction term is the same. 

With the aid of these independent counterterms, we finally arrive at the following 2PI 
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(a) 



(fa) 



FIG. 2. 2PI diagrams at the two-loop order: (a) the double bubble diagram and (b) the setting 
sun diagram. The double bubble diagrams has the contribution of the s-channel. On the other 
hand, the setting sun diagram has the contribution of the t- and n-channels. 
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FIG. 3. The counterterm diagrams of (a) the tadpole consisting of the contraction of the two point 
interactions, and (b) the double bubble consisting of the contraction of the four point interaction. 

effective action at the 2PI Hartree approximation in the symmetric phase 



r[v\G, 



d X 



-vUD + mi)v 



2\^.i , '^4 / .i .iN2 
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\v V 
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{U + ml + -\^{Y. 



^'^') + ^Afi;V)G',,(x,x') 



(9) 



d'x k E ^"(^' ^) E ^^•^■(^' ^) + 2^^ E ^^^(^' ^) 



where we have introduce the independent couphng counterterm 5A4 for S[v\ as in the case 
of 0^ theory and m^ = m^ + 5mQ. 

B. Renormalization in the broken phase 



In this paper, we are interested in the broken phase and we renormalize the ultraviolet 
divergence in the broken phase. Of course, this renormalization is consistent with that in 





5\i + 2S\f 



5X4 




SX^ + 2S\B 





(JV - 1)5A^ + 2SXB 



FIG. 4. The counterterm diagrams at the Hartree truncation in the broken phase. The double Une 
represents the a full propagator and the single line represents the vr* full propagator. 



the symmetric phase. 

We now break the phase by specifying the vacuum expectation values as v^ = v, v 
0, z 7^ 1 and rename the fields as v?^ = cr, v?* = tt*, i 7^ 1. In these variables, the 2PI effective 
action is reexpressed as 



T[v,G^,G^] 



d X 



-v{n + mo)v + —V 



+ ^lndet[G;^] + ^(iV-l)lndet[G;i] 



'-NTr[l] - iTr[(n + mg + 1(A^ + 2X^^)0^] 

i(iV - l)Tr[(n + ml + ^A^^)G.] - ^ |rf^x(A^ + 2X^)Gl 

{N-1)X^ + 2X^' 



^(N - 1) fd'xXtG.G^ - ^{N - 1 



d X 



Gl 



(10) 



where we have identified the various species of the fields vr* as tt. The counterterm dia- 
grams contributing to this expression are shown in Fig. HI In the above expression, the 
different counterterms derived by the different contractions effectively appear as different 
counterterms for the a and vr* fields. 

Equations of motion for f , Gg-, G^^ are derived by the derivatives of the 2PI effective 
action with respect to these variables. They are written as follows: 



D + m^ + ^X,v^ + ^(^^ + 2^^)G. + ^(^ - ^)^2G. I v{x) 



0, 



fill 



U+ml + ^{X^+2X^y + ^{X^+2X^)G„ + ^{N-l)X^Gj\G,{x,x') = t5{x-x'), (12) 

U + ml + 1a^2 + l_x^G^ + ^ [(^ - 1)^^ + 2Af ] G^] G^{x, x') = t5{x - x'). (13) 

We first assume that the equations of motion are renormalized by the MS scheme, that 
is, only the divergent terms are removed by the counterterms, and derive the equations 
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for the counterterms to satisfy. Then, we exphcitly construct the counterterms by solving 
them. From the above equations of motion, we see that one needs 6Xq ' = 6X2 ' and 
6X4 = 6X2 + 26X2 for a consistent renormahzation. Strictly speaking, the last equality can 
be derived from the renormahzation of the four point vertex function of the a field. This 
can be seen in Appendix A. 

From the equations of motion for the propagators, we define the mass of the a and vr* 
fields as 



ml = ml + i(A^ + 2Xiy + j^{X^ + 2X^)G^ + j^{N - l)A^a, 



m. 



771 



+ ^A^2 + 1 A^G. + UiN- 1) A^ + 2X^] G^, 



N ' N "^ N 

and assume that these expressions are renormalized as follows 



(14) 
(15) 



ml = m' + ^Xv' + |at; + l(iV - 1)AT;, 



(16) 



ml = m' + ^Xv' + Iat; + ^{N + 1)XT^, 

where we have used the divergence structure for tadpole correction, G'o-,7r = m^^Ta 
Then, the counterterms have to satisfy the following equations. 



(17) 






67710 + j^{6Xt + 26X^y + ^(3A + 6Xt + SX^)mlT, + 1(<5A^ + 26X^)T^ 

+ l(Ar - 1)(A + 6X^)77ilT, + l(Ar - 1)6X^T^ - 



(18) 



+ ^[{N + 1)A + {N- 1)6X^ + 2dX^]7nlT, + ^[{N - 1)6X^ + 26X^]T^ = 0. 



(19) 



We can explicitly construct the counterterms for the MS scheme by solving these equations. 
For this, a crucial step is to use the renormalized expressions for the mass. Collecting the 
terms proportional to f^, Tp and Tp, we obtain 



67nn 



^{N + 2)77i^XTd + 6X^m^Td + pX^77i^Td 



+ v^ — 

N 



+ TZ- 



N 



^{N + 8)X^T, + 6X^ + 26X^ + ^(^ + 2)a^AT, + ^^0""^^^ 
^{N + 8)X^Td + SX^ + 2SX^ + ^(^ + 2)a^Arrf + ^SX^XT, 



+ T^-{N-1) 



-{N + 4)X'T, 



S^t + ^(^ + 2)5A^T, + ^5A^AT, 



(20) 



5mo + ^iN + 2)m^\Td + SX^m'^Td + ^SX^m^Ta 



+ v' 



+ t; 



+ TZ 



N 
1 

N 
1 



i(iV + A)X'T, + 5Xt + ^{N + 2)5X^XT, + |a^AT, 



(21) 



N 



1 

iV' 



(iV + 3)A2t, + (iV - 1)5A^ + 25A^ + ±.(N-l){N + 2)5X^XTd 



+ -(N+1)6X^XT, 



0. 



Here, we assume the divergent terms proportional to v"^ and Tp as sub-divergence which 
comes from the divergence of sub-diagrams and the residual divergent terms as overall- 
divergence. 

We then assume that the terms proportional to v"^, Tp and Tp independently vanish. The 
vertex counterterms 6X can be determined by the cancellation of the sub-divergences, and 
the mass counterterm 6mQ is determined by the cancellation of the overall-divergence. In 
fact, we obtain the explicit expressions for 6X from the equations of Tp and Tp 



SX^ 



-X^T 



N 

oo 



2 1-1 



n=l 



^A^;^ 



l + l(iV + 2)AT;"' 



A^ 



AT, 



(A^ + 4)A + 25A^ 



(22) 



(23) 



The mass counterterm 5mo is constructed from the cancellation of the overall-divergence as 
Smn = — m^Tij 



1-^N + 2)X + 5X^ + ^5X^ 



m^ — (N + 2)XTd 



N 



l + -(iV + 2)AT, 



(24) 



-m 



E 

n=l 



N 



(iV + 2)AT, 



We have renormalized all of the equations of motion. The effective action is renormalized 
by these counterterms. 

When we work on curved space, we have to add the following coupling term to background 
geometry in the action for a consistent renormalization. 



S[4>i] = - jd^xyf-g 



]^ct>\U + m^ + iR)4>' + ^{ 



kiA.i\'2 



(25) 
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where R is the Ricci scalar and ^ is a numerical factor of the conformal parameter. The 
alteration in the equations of motion is the coincident propagator. In a general curved 
background, it is found by the heat kernel technique that a coincident propagator, which is 
regularized by the dimensional regularization scheme in four dimensions, has the following 
general form [l5|. 



G{x,x)= [m' + {i-QR]Td + Tp, 



(26) 



where ^c = {d — 2)/4((i — 1) and Tp is an arbitrary tadpole correction. This alteration 
gives rise to further overall ultraviolet divergence (^ — ^^^RTd in the equations of motion. 
However, the sub- divergence terms of w^, Tp and Tp undergo no alterations. That is, the 
sub-divergence terms are renormalized by the same counterterms 5\ as fiat space. The 
additional overall-divergences are renormalized by the counterterm of the newly introduced 
conformal parameter 5S,o- It has the following form similar to the mass counterterm. 



5io={i-QT, 



i_(Ar + 2)A + a^ + |5Ao 



B 



-{i-Q-{N + 2)XTa 



-\ -1 



l + -(iV + 2)AT, 



(27) 



oo . V n 

n=l ^ ^ 



The effective action is also renormalized by these counterterms and the redefinition of the 
coupling constants in the gravitational action with curvature squared terms, in a similar 
way to 0^ theory 13 1. 



III. NONPERTURBATIVE INFRARED EFFECTS IN DE SITTER SPACE 



Throughout this paper, we work with the standard in-out formalism. In a certain curved 
background, like de Sitter space, the metric has a time dependence, and its nonequilibrium 



ism 



nature may appear. In such a situation, it is known that the standard in-out forma 
is not sufficient, and it is more appropriate to take the Schwinger-Keldysh formalism 16 |. 
However, we omit the closed-time path index for the Schwinger-Keldysh formalism, since for 
our approximation order, at the Hartree truncation level of the 2PI effective action, these 
in-in and in-out formalisms give the same results. 



11 



We use the coordinate system for de Sitter space in terms of comoving spatial coordinates 
X and conformal time — oo < rj < in which the metric takes the form 

= a{T]f{drf-dy.^), 
where 0(77) = —l/Hrj is a scale factor and H is a Hubble parameter constant. 

From the renormalization in Sec. II, the renormalized effective action is given by 



(28) 



T[v,M^,M^] 



+ 



3 

In 



2 ^ ' 



A^4 _ 1 / ^^2y. 

AN 2 



V-^F 



2^ 



dMlT^p 



1 

2iV" 



XTf + -^(iV - 1)AT^T; + -^(iV - l)(iV + l)\Tp 



1 

AN' 



•■k2 



(29) 



Furthermore, the renormalized mean field equation and the self-consistent mass equations 
are respectively given by 

w = 0, (30) 



Ml - -Xv[ 



Mi = M^ 



Mi = M' 



^XT^ + 1(N-1)\T^, 



^Xv^ + ^Xn + ^{N + l)XT^,, 



(31) 
(32) 



M' 



where we have assumed that 1; is a constant due to the spacetime symmetry of de Sitter 
space. From these equations, we can investigate the possibility of spontaneous symmetry 
breaking. To do so, we first assume that v has nonzero solutions and M^ — 2Xv'^/N = 0. 
Using this equality, we eliminate the v dependences in the self-consistent mass equations. 
Then, the self- consistent mass equations are given by 

N ^ N 
2 2 

The possibility of spontaneous symmetry breaking can be investigated from whether these 
equations have solutions. Now, the tadpole correction Tp is expressed as the digamma func- 
tion (see Appendix B). An analytic search for solutions requires an approximate expansion 
for the tadpole corrections Tp and Tp. In doing so, we first assume that the Goldstone 
bosons are light compared to the Hubble parameter H. Then, the tadpole correction Tp is 
expanded as 



-2M^ 



Mi 



:(iV-l)AT-, 

:AT? 



F- 



(33) 
(34) 



T^ 



167r2 



M2 



+ o{i^ 



(35) 
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Next, the expansion for Tp is possible in two ways. One is to assume that M^- is heavy 
compared to H. Then, the tadpole correction Tp is expanded as 



n 



H 



2 



^ 167r2 



3 ^^Im^ 



(36) 



where we have set the renormalization scale as — 1 + 7 + log(^^) = 0. These expansions 
enable us to solve the self-consistent mass equations as an algebraic equation. Taking the 
lowest order of the tadpole expansion, the solutions are respectively given by 



M: 



2 



..2 ^H' 1/ \H^- \2 3\H' 

These solutions indicate that our approximate expansions, M^/H"^ ^ 1 and M^/H"^ ^ 1, 
are justified in the parameter region A <^ 1 and H^/\M^\ <^ 1, i.e. weak coupling and weak 
curvature. The solutions reproduce the masses with spontaneous symmetry breaking in flat 
space in the limit H -^ 0. Thus, these solutions correspond to the masses at the minimum 
of the effective potential. The tadpole correction Tp can also be approximated by the small 
mass expansion. In this case, the self-consistent masses are given as the solutions of an 
algebraic equation of the fourth degree at the lowest order of the small mass expansion. 
This approximation is discussed in the earlier study (l4l |. 

The above analysis means that spontaneous symmetry breaking is possible in de Sitter 
space. Also, the mass system has two solutions at some Hubble parameter. This means 
that the effective potential has two extrema, and shows the signal of the first order phase 
transition. From Eq. ( l38l) . the mass of the Goldstone bosons cannot be zero. That is, the 
infrared divergence in the propagator is self-regulated. 

These analytic results are reliable only in the limited range of the parameters. In seeing 
the behavior of the solutions in all parameter regions, we have to resort to numerical cal- 
culation. Numerical solutions for the mass system of Eqs. (133|1 and ( 134|) are shown in Fig. 
5, where the renormalization scale is set as —1 + 7 + log ( '^^^2 ) = 0. The solution in the 
symmetric phase t' = is also plotted. The figure shows the critical value of the Hubble 
parameter He at which the broken phase solutions become to disappear. This indicates the 
phase transition at He- Moreover, the mass system has two solutions at some value of H. 
This is a signal of the first order phase transition. The mass of the Goldstone bosons cannot 
be zero. These results are consistent with those obtained by the analytical method. 
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FIG. 5. The solution of the mass gap system as a function of a Hubble parameter H in the case 
of M^ < 0. A^ = 6, A = 0.8 and all values are in the units of |M|. The upper and lower curves 
represent the mass with v on the maximum and minimum respectively of the effective potential as 
shown in Fig. [6l The straight line represents the solution in the symmetric phase. 



The effective mass term of the Goldstone modes superficially appears to be a violation 
of Goldstone's theorem. However, we can attribute the mass generation to the lack of 
time-translational invariance in de Sitter space [l7|. In such a circumstance, spontaneous 
symmetry breaking does not necessarily imply the massless Goldstone bosons. Thus, the 
effective mass generation is compatible with symmetry breaking and Goldstone modes. 



Here, the following comment should be made. The solutions of the mass system in 
Fig. 5 are similar to the same analysis at finite temperature 18|. In fact, it was found 



in various contexts that de Sitter space has thermal aspects with the temperature given 
by Th = i//27r pJiJ. This is due to the existence of the cosmological horizon of de Sitter 
space. However, Fig. 5 also shows the difference between finite temperature. A difference 
is that the lower curves of the mass solutions always start at iJ = 0. In the terminology 
of condensed matter physics, the lower spinodal point iJ^i always exists at zero value. On 
the other hand, in finite temperature field theory, the lower spinodal point exists at nonzero 
values of temperature. 



The nature of the phase transition indicated by the analysis of the mass system can be 
consistently confirmed by plotting the shape of the effective potentials. From the renormal- 
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FIG. 6. The effective potentials as a function of v near the critical temperature for A^ = 6 and 
A = 0.8 all in the units of \M\. The different lines show the potentials with different values of H. 
Normalization is taken so as to match the origin for different values of H. 

ized effective action Eq. ( !29l) . the effective action is given by 

" W^^^' ~^^^' ^W^fT^f - ^{N - l)(iV + l)\Tf. 
Numerical calculation of the effective potentials near the phase transition is shown in Fig. 
6, where the renormalization scale is set as — 1 + 7 + log( ~^^ ) = (see Appendix C). In 
fact, the first order phase transition takes place at the value of Hf. which is indicated by 
solving the mass system. 

In the earlier study of the large- A^ expansion, the strong infrared effects prevent the pos- 
sible existence of spontaneous symmetry breaking |5| . In other words, infrared divergences 
of the Goldstone modes are circumvented by prohibiting the spontaneously broken phase. 
This fact is reminiscent of the absence of the spontaneous breaking of continuous symmetry 
in fiat two dimensions. In contrast, our result demonstrates that infrared divergences of 
the Goldstone modes are circumvented by generating the effective mass term of the Gold- 
stone modes rather than the absence of symmetry breaking. These results appear to be a 
contradiction. However, the conflict is reconciled by observing the behavior as A^ increases 
in our analysis. The behavior of the mass solution M^ as A^ increases is shown in Fig. 7. 
The figure predicts that the critical value of H decreases as A^ increases, eventually He goes 
to zero as A^ goes to infinity. In fact, if we assume that the small mass expansions of the 
tadpole corrections Tp and Tp are valid near the phase transition, the critical value of H 

15 




0.0 0.2 0.4 0.6 0.8 
H 

FIG. 7. The solution of the mass gap system for M^ as a function of a Hubble parameter H in 
the broken phase for various number of A^. A = 0.8 and all values are in the units of \M\. These 
solutions indicate that the critical temperature He decreases as N increases. It is anticipated that 
the critical temperature goes to zero in the limit A^ — )■ oo. 

can be obtained analytically at large value of N . At the leading order of 1/A^, it is given by 

/ 1 V* 

The validity of this expression is confirmed numerically. From the expression, we see that 
He becomes zero as A^ goes to infinity. This behavior is different from that at finite temper- 
ature, there the Goldstone bosons become massless as A^ goes to infinity. When the critical 
curvature is zero, spontaneous symmetry breaking is prohibited and the broken phase in 
flat space undergoes a discontinuous transition. Thus, it is concluded that the large- A^ ex- 
pansion only tells the nature of the theory as A^ goes to infinity, predicting the zero critical 
curvature. 

IV. CONCLUSION 

We investigate nonperturbative infrared effects for the 0{N) linear sigma model in de Sit- 
ter space. Ultraviolet divergences are regularized by the dimensional regularization scheme. 
In all calculations, we work with coordinate space representation. 

As a results of our analysis, spontaneous symmetry breaking is possible for any finite 
number of A^. In the broken phase, the possible infrared divergence is circumvented in a 
manner of generating the mass term for Goldstone modes. The broken symmetry is restored 
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by the first order phase transition as a function of the Hubble parameter. These results have 
a similarity to finite temperature field theory with the temperature given by Th = H/2'jt. 
A difference with finite temperature field theory is that the critical temperature has the A^ 
dependence of H^. ~ N~'^^^, predicting a zero value of the critical temperature as A^ goes to 
infinity. From the result, we conclude that the large- A^ expansion only tells the behavior of 
the 0{N) model as A^ goes to infinity, i.e. the absence of spontaneous symmetry breaking 
and the discontinuous transition. Thus, our results are compatible with those obtained by 
the large- A^ expansion. 

The effective mass term of the Goldstone modes superficially appears to be a violation 
of Goldstone's theorem. However, we attribute the mass generation to the lack of time- 
translational invariance in de Sitter space. In such a circumstance, spontaneous symmetry 
breaking does not necessarily imply the massless Goldstone bosons [17|. Thus, the mass 
generation is compatible with symmetry breaking and Goldstone modes. The effective mass 
term for Goldstone bosons is also predicted in the same analysis at finite temperature, and 



this is generally considered as an artifact of the Hartree approximation [12|, |20[ . However, 
we believe that the mass of Goldstone bosons at finite temperature is also caused by the 
lack of time-translational invariance at the formalism in treating the finite temperature, like 



the imaginary time formalism 211 ]. 



All our results are based on the mean-field approximation. It is necessary to include some 
contributions of the non-local diagrams, such as the setting-sun type, for a more rigorous 
study. For example, in a similar study at finite temperature, it is found that the order of 



the phase transition 



may 



_3e changed from the first to the second by the contribution of 
non-local diagrams 22, l23|]. This is because in our analysis, it is shown that the results 
in de Sitter space does not necessarily share the same nature with finite temperature, like 
its behavior as A^ varies. Furthermore, the recent study of the 0(2) model by the Wigner- 
Weisskopf nonperturbative method also predicts the mass generation and first-order phase 



transition 17l |. 



The mass generation of Goldstone bosons is stated in another way that the self-regulation 

n 

mechanism for IR divergence mentioned in Ref. [6| has taken place in the more complex 
model. We anticipate that the self-regulation mechanism will take place in quantum gravity 
on a de Sitter background due to field self-interactions. 
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Appendix A: Four-point vertex function of the a field in 0{N) model 

In this Appendix, we confirm that the renormahzation condition 6X4^ = 6X2 + 2(5Af 
which comes from the renormahzation of the equations of motion is consistent with the 
renormahzation of the four point function of the a field. The four point vertex function of 
the a field is derived as a fourth derivative of the IPI effective action as the general form: 

5^T S^T 



5'ripi 



S'Gl, 



6V16V2SV36V4 



6V16V2SV36V4 



Svi6v2SG\ SV3SV4 



+ (6 perm.) 



+ 



''^ "^ + (Sperm.), 



(Al) 



^G\bG^^^ bvxbv2 bv^bv4 



where we have used a short-hand notation v\ = v{xi) and G12 = G{xi,X2), "perm." denotes 
the possible permutations of the indices 1, 2, 3 and 4, and the repeated indices i and j run 
from a to TT. Here, we use the relation 



6^ 



rm r ITT /^TT— 1 /^TT — 1 



SG:,SGi, 



'-S'''6^''g:~'gz' + ^(iv - i)6'-6^''g:;'gi. 



S^T, 



SGl,SGi, 



(A2) 



and 



S^G 



ab 






(A3) 



5vi5v2 

where the index {€) denotes the subsidiary index following to the repeated index i. Then, 
the four-point vertex function is transformed to 



^(4) 
IPI 



5^T 



53r 



,5'^G 



i-l 



5vi5v26v^5v4 5vi5v26G 



ab 



6V2SV4"^'' 



/^{i) ^ ^^ cd /^ 



(6 perm. 



i S G^j- 



f, 



5'^G 



0--1 
dh 



"^f-TT— 1 



GfJ 



2 5vi5v2 ^^ SV3SV4 ^^ 



5^T, 



,<52G 



i-l 



SGUSGi, 



g^I—^gLg^^} 



5^G- 






(A4) 






5vi6v2 '^^ SV35V4 



G'^> + {Sperm. 



We now rewrite the last term by the Bethe-Salpeter equation for S'^G ^/Sv"^. The Bethe- 
Salpeter equation for the a field is derived as follows 



(^ Ggb 

6vi6v2 



6' 



6vi6v2 
■■-2i 



Go^ - 1i 



bV2 

sg:. 



6^T, 



5^G 



cd 



-2i 



SG:,6v^6v2 6Gl,5Gl,5v,5v2 



5^V2 



5G:,5v,5v2 6G:,SG 






(A5) 
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In the same way, the Bethe-Salpeter equation for the vr field is 



6' 



6vi6v2 6vi6v2 
2i 



Gq 



2i 6T^ 



N-1 



N- 



l^G-,, 



^^r, 



SG:,M6v2 6G:JG 



cd 



Mi) e/ ^{i) 



(A6) 



Using the Bethe-Salpeter equation, we obtain 



.(4) 
IPI 



5^V 



5vi5v25v'i5v4 



5^T 



-GS?S^g2 + (Sperm. 



(AT) 



The last term can also be rewritten by the Bethe-Salpeter-like equation. To see this, we 
write down the various four-point kernels. 



5^T 



SvSv6Gc 

6^r 



-^(A^ + 2Af), 



SvSv5Gn 

^2 



S^T, 



1 

N 
1 



-^(A^-l)A^, 



5Gfj5Gfj 

^2 



^^r, 



5GrrSGT, 

^2 



6^r, 



2N 
1 

'2iV 
1 



(A^ 



2Ao^), 



(Ar-l)A^, 
(iV_2)[(iV-l)A^ + 2A^ 



(A8) 

(A9) 

(AlO) 

(All) 

(A12) 



SGJG^ 2N' 

where we have omitted the spacetime indices. We find that the last term is transformed by 
the Bethe-Salpeter equation for b'^G"^ jbv^ . Finally, the four point vertex function of the a 
field is expressed as 



.(4) 
IPI 



5^r 



bv\bv2bv-ibv4 

6 , 6 

N" ' N 



6^r, 



1 6'Gl,' 



SGi25v3Svi 2i SV3SV4 



(6 perm. 



ttA + ^(-(5A4 + 5A^ 



2S\i) + - 



^-|- + -A)+(6perm.) 



(A13) 



6\o + 26X2 is required for the 



From this expression, we see that the condition 6X4 — u/\2 

renormalization of the four point vertex function of the a field. This condition is consistent 

with the renormalization of the equations of motion. 



Appendix B: Coincident propagator 

In this Appendix, we calculate the coincident propagator in de Sitter space to investi- 
gate the divergence structure of the tadpole diagram. In d dimensional de Sitter space. 
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the propagator for a free scalar field with mass m, conformal factor ^ is expressed by the 
hypergeometric function 24 1 



G{x,x') 



jjd~2 r(^ + zy)r( 



d-i 



V] 



(4vr)rf/2 



2Fi[^ + i.,^-i.,f;l + |], (Bl) 



r2\i/2 



where v = {[{d — l)/2]^ — (m^ + ^R)/H'^} , R = d{d — l)H'^ is the Ricci scalar curvature 
and |/(x, x') = [(77 — rj'Y — |x — x'p] /rjr]' is the de Sitter invariant length. In the coincident 
limit, y = 0, the formula of the hypergeometric function, 2Fi(a, 6, c; 1) = T{c)T{c — a — 
b)/ [r(c — a)r(c — b)~\ , leads to 



G{x,x) 



H' 
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(4vr)d/2 
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-rn 



./(¥ 
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^) 


^^ r(| 
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^)r(| - 


^) 




f)r(x,x 


)• 









(B2) 



The first gamma function has an ultraviolet divergent pole. The residual gamma function, 
r(x,x) determines a coefficient of the ultraviolet divergent pole. 

We restrict our attention to four dimensional spacetime with a regularization parameter 
e = 4 — d. In this case, we can transform the expression r(a;, x) as follows: 



T{x,x) 



r(i + 



d-3 



u)T{l + ^-u) 



va+uni-u 



d-3 



2 
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2 y V 2 y H^ 

r(i + z.) [1 + ^(1 + z.)(-f ) + 0(e2)]r(i - z.) [1 + ^(i - z.)(-|) 



o(. 



m^ + iR 



-id-T 



^A{^' + {^-m 



r(i + z.)r(i - z.) 






\i.R 



(B3) 



where ip{x) is the digamma function and ^c = {d — 2)/4((i — 1). Therefore, the coincident 
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propagator is generally given by 



G{x,x) 



167r2 



^ { (^' + (e - QR) [i - (I) (v^(i + ^) + ^(i - ^)) + 0(e^) 



rCci? , 



167r2 
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(B4) 

where 7 is the Euler-Mascheroni constant. The renormalization scale /i is introduced in 
advance. Strictly speaking, it is introduced when the coupling constant of the field interac- 
tion is made dimensionless in the dimensional regularization scheme. The above ultraviolet 
divergence structure appears in any curved background. In fact, it is revealed by the heat 
kernel technique that the coincident propagator in an arbitrary background has the form 



G{x, x) = [m^ + (e - ^c)R]Td + Tf, 



(B5) 



where T^ = —2/lQ7r'^e and Tp is a general finite tadpole correction. That is, the finite 
tadpole correction Tp in de Sitter space is given by 



Tf = Y^h'+(e-ec)i?] (v'(|+^)+^(|-^)-i+7+iog^+iog^ 



MM 1 

Note that the above expressions correctly reproduce the fiat space results in the limit iJ — ?> 0. 
In the fiat space limit, canceling the H dependent terms, the coincident propagator is given 
by 



G{x,x] 



flat 



771 
16^ 



1 + 7 + log 



m 



ilTfl^' 



(B7) 



In performing various calculations, we have to evaluate the tadpole correction Tp. The 
digamma function in the tadpole correction can be approximately expanded regarding the 
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magnitude of the mass term M^. First, in the small mass case M"^ /H^ ^ 1, Ti? is expanded 

as 



T, 



16^ 



6H^ 20 



H^ 



M2 



_|,4.-2(-l+-, + log£^)H-0,fJ 



(B8) 



On the other hand, in the large mass case M'^/H'^ ^ 1, Tp is given by 



Tp 



M 
16^2 






AH^ 



47r/iV 3M2 
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Appendix C: One-loop effective action 



In this Appendix, we calculate the finite part of the one-loop effective action. This is 
expressed as the mass integral of the tadpole correction Tp, 



nXop = -\j d'x^^ jdM'Tp{M^). 



(CI) 



Using the result in Appendix B, we obtain 



dM'^TpiM'^ 
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where 0{e) terms are omitted. Also, renormalization scale fi is introduced in advance. 
The integral of the first term can be performed by changing the integral variable M^ to 

1/ = (9/4 - MVi72)i/2, as follows 



.1 



=z/(--z.2)[logr(l + z/)-logr(i-z/)]-(i-3z.2)[V^(-2)(i + z/)+V^(-2)(i-z/)] 



(C3) 
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Finally, the finite part of the one-loop effective action is given by 



dM^Tp 
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iTT" 
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The above expression has the correct fiat space limit as H goes to zero: 
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